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DIFFUSION WITH CHEMICAL REACTION: AN ATTEMPT TO EXPLAIN 
NUMBER DENSITY ANOMALIES IN EXPERIMENTS 
INVOLVING ALKALI VAPOR 

By Walter L. Snow 
Langley Research Center 

SUMMARY 

This report theoretically examines the mutual diffusion of two reacting gases which 
takes place in a bath of inert gas atoms. Solutions are obtained between concentric 
spheres, each sphere acting as a source for one of the reactants. 

The model is an attempt to explain recent anomalies in alkali densities which have 
been observed in optical absorption experiments. Severe spatial gradients of alkali den- 
sity were noted when uniform distributions were expected. This paper suggests that the 
observations might be explained in terms of the alkali reacting with outgassing and desorp 
tion products from the walls of the vessel. 

Coupled equations for spatial density distributions of the reacting components are 
provided; the coupled equations can be solved by iteration and depend only on diffusion 
coefficients, second-order rate constants for the reaction, and geometric dimensions of 
the apparatus. Formulas are provided to estimate the rate constant when it is otherwise 
unavailable. Results should be applicable to the treatment of other chemically reacting 
species. 


INTRODUCTION 

Recent measurements (ref. 1) of the alkali number density distribution over reser- 
voirs of the molten metal have uncovered very severe spatial gradients when uniform 
distribution is expected from cursory kinetic theory considerations. Typically, the 
reservoir of alkali is placed in an isothermal vessel containing several hundred torr 
(1 torr = 133 N/m2) of inert gas. Mean free paths and times for alkali-inert gas 
encounters are, respectively, on the order of 10" 5 cm and io-9 sec at the temperatures 
considered. In seconds the alkali density is expected to be uniform throughout the con- 
tainer and to have a value comparable with that predicted by using vapor pressure data. 
Observed variations of over a decade change per cm are not uncommon. Figure 1 
describes typical experimentally observed alkali number density distributions. 



The vessels used in the experimental measurements of reference 1 were prepared 
with moderate care. Typically, after rough pumping for several hours with mild heating 
to accelerate outgassing, the alkali was added in a nitorgen-flushed glove box. The ves- 
sels were flushed with dry inert gas and subsequently closed, pumped, and backfilled to 
several hundred torr of inert gas. However, some residual desorption and outgassing is 
still to be expected under these conditions. 

This paper theoretically investigates whether residual desorption and outgassing 
levels in the early phases of experimentation (usually several days) could account for the 
observed spatial gradients. The calculational model is used to illustrate the types of 
effects observed in reference 1. 

The geometry of the experimental setup is dictated by considerations for measuring 
optical transmission. The source of alkali is preferably kept on the periphery of the con- 
tainer and outside the light path. In the present study the theoretical model is treated in 
a simpler, symmetric configuration to facilitate solving the related partial differential 
equations with boundary conditions. It is felt that the similarities will overshadow the 
geometric differences enough to make the description applicable to experiments in 
reference 1. 

The idea for treating this problem as diffusion with absorption (through reaction), 
was considered as an analog to the thermal diffusion of neutrons in a reactor with modera- 
tor. (See ref. 2.) Girifalco and Behrendt (ref. 3) also treat the diffusion equation with 
source and sink terms in spherical coordinates. They employ eigenfunction expansion 
techniques and have source terms distributed throughout, rather than being localized on 
the boundaries as in this case; therefore, their results are not applicable here. 

The consequences of the effect described in this paper can be drastic in those situa- 
tions where the gas phase is calculated from vapor pressure data on the liquid or solid 
phase of a reservoir whenever a reactive constituent is present even in trace amounts. 

In closed systems where there is no liquid reservoir of material present, the reaction 
would, of course, quickly deplete the reacting gas phase constituent. In other situations 
when gradients exist and average densities are calculated on the premise of constant den- 
sity, measured density values would appear to be lower than those predicted from vapor 
pressure equations. Reference to both of these effects (sometimes under the guise of 
nonrepeatability) is frequent in experimental papers concerned with absorption vessels 
and discharge tubes. The high reactivity of alkali makes the effect easier to observe 
than usual. 
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SYMBOLS 


A species type, mnemonic for alkali 

B species type, mnemonic for buffer or bath 

C species type, mnemonic for contaminant 

Cj constant (i = 1, 2, . . ., 6, cm _ 2 

D diffusion constant, cm2- sec - 1 

Djj generalized diffusion coefficient, cm2-sec“l 

binary diffusion coefficient, cm2-sec _ l 
AE activation energy, erg 

f( ) auxiliary function defined prior to equation (9), cm- 2 

G Laplace transform of g, cm - 2-sec 

g( ) auxiliary function, = tt ], cm "2 

g parameter modifying annihilation constant, dimensionless 

h Planck constant, 6.63 x 10-27 erg-sec 

h(r) auxiliary function defined prior to equation (12), cm-2 

I moment of inertia, g-cm2 

J flux, number of particles per unit area per unit time across a surface whose 

normal is specified by the vector, cm - 2- sec -1 

k Boltzmann constant, 1.38 x 10-23 j-K'* (l.38 x 10“ erg-K“l) 

Laplace transform operator 
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N(T) 


n 

© 

P 

Q 

q 

R 

r 


s 

T 

t 


X' 


z 


inverse Laplace transform operator 
mass of species i, g 

reduced mass of species i and j, m^m^m^ + m.^, grams 

number density obtained from vapor pressure (or outgassing) for temperature 
T, cm~3 

number density, cm -3 

order of magnitude 

pressure, torr (1 torr = 133 N/m2) 

collision cross section, cm2 

time constant for outgassing rate, sec 

radius of container, cm 

radius coordinate, cm 

radius of source bead, cm 

Laplace transform partner for t, sec-1 

temperature, K 

time, sec 

time of impulse perturbation in source strength, sec 
activated complex concentration, cm -3 
partition function for activated complex, cm-3 
partition function for single mode, dimensionless or cm-1 
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0 


attenuation parameter, cm _ l 
y attenuation parameter, cm - l 

A perturbation strength parameter, dimensionless 

6( ) Dirac delta function 

N. 

e ratio of vibrational to rotational partition functions per mode, dimensionless 

£ attenuation parameter, cm - 1 

7 } number density excess, cm - ^ 

k annihilation constant, sec'l 

v frequency, sec _ l 

t = t - t Q , sec 

Superscripts: 

d dissociation 

e electronic 

r rotational 

t translational 

v vibrational 

°o steady state 

Subscripts: 

a alkali 

b buffer 
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c 


contaminant 


i,j,k labels for arbitrary species or constants 

* activated complex less one degree of vibrational freedom 

THEORY 

The model for this problem is sketched in figure 2. A small bead of alkali metal 
having radius r Q is situated in the center of a spherical cavity whose radius R is 
usually much greater than r Q . The concentration of gas phase alkali is of interest and 
is designated n a (r,t). Instead of treating the heterogeneous phase problem from first 
principles, it will be assumed that the liquid phase (temperatures of interest are above 
the melting point of the alkali) is in equilibrium with its vapor in a thin shell surrounding 
the bead. The concentration of gas phase atoms at r Q is taken from equilibrium vapor 
pressure data; thus, the difficult interface problem is avoided. 

Atoms of type A (alkali) diffuse outward into a bath containing predominantly spe- 
cies B (buffer gas). Also present are contaminant atoms designated C which desorb 
and outgas from the walls at R. Atoms B are inert and react with neither A nor 
C, but encounters of A and C can result in the formation of AC, which effectively 
removes A (alkali) as an atomic constituent. 

An alkali atom leaves the vicinity of the bead and moves at several thousand meters 
per second. The alkali atom is constantly battered by B atoms; this condition makes 
A atom's path very tortuous and improves its chances for an encounter with C before 
it can reach the outer wall. If only A and C were present at their low partial pres- 
sures, encounters would be very unlikely and, in fact, the mean free paths would be gov- 
erned more by the container dimensions than by atomic encounters (Knudsen regime). 

The mutual diffusion of A and C will result in gradients of n a because of chemical 
reaction. The presence of B is acknowledged through use of a generalized diffusion 
coefficient to be defined later. 

A comprehensive and readable discussion of diffusion is to be found in reference 4. 
For the homogeneous phase problem the net flux of particles is governed by Fick's 
empirical law 

J = -D(r ) Vn(r,t) (1) 

where J is the net number of molecules crossing unit area per unit time in the direction 
specified by J, n is the number density, and D is the diffusion coefficient. In a given 
volume element the only mechanism for n to change is by chemical reaction; the reaction 
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rate is denoted by an annihilation constant k, which usually is expected to depend on the 
concentration of reactant present. (See the following equations and the appendix.) 
Mathematically, 

( 2 ) 

dt 

or 

^^+9* J = -xn(r,t) (3) 

By using Fick's law, 

8n l+*^ ~ VD(F) • Vn(F,t) - D(F)V 2 n(F,t) = -/cn(r,t) 
at 

Rearranging this equation yields 

8n ^’ t -' = D(r )V 2 n(r,t) + VD(r) • Vn(r,t) - xn(r,t) (4) 

In practice, the diffusion coefficient is governed by the bath atoms and its gradient is 
negligible. Henceforth vB(r ) is assumed to be zero. An explicit expression for the 
annihilation constant is derived in the appendix. It is reasonable to suppose that it will 
depend functionally on the concentration of reacting species. Solutions to the following 
coupled equations are sought: 

ana 9 ( t r,t) = Dac^^a^ ' K a( n c) n a (F ^ < 5a ) 

9nC ^ r - * D ac (r )v 2 n c (r,t) - * c (n a )n c (f,t) (5b) 

The experiments in reference 1 took place over several hours or days. Mostly, it 
is the steady-state solution which is sought, and these concentrations are identified with 
a superscript n(r,t=<=o) s n°°(r). An approximate method of studying the damping of per- 
turbations to the steady-state values is discussed subsequently with possible application 
to empirically estimating «. 
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Steady -State Solutions 


Possible anisotropy in the problem (for example, gravitational effects) is 
neglected and spherical symmetry in alkali-inert atom collisions is assumed. Then, 
n(F,t) = n(|F|,t) = n(r,t) and the Laplacian operator simplifies to 


r;2 _ 1 d r 2 d 

r 2 dr dr 


( 6 ) 


Expression for n^Cr).- The boundary conditions for the problem will be 


n a* 9 ( r o) * N a(T) 

and 


(7a) 


vi^(R) = 0 (7b) 

The first equation follows from restricting the number density to be the equilibrium value 
at the interface to sidestep the multiphase aspects. The second equation follows from 
Fick's law and the requirement that there be no net flow of alkali across the container 
boundary. Equation (5a) reduces to 


D ac ^ r2 SF "a" (r > - l ‘a( n c) n a"M = 0 


( 8 ) 


The equation for the related function f = rn °°(r) is even simpler, for example, 

a 


&-V 2 f = 0 
dr^ 


(9) 


with r 2 = K a (n c y Dac . 

Equation (9) is nonlinear and is coupled to equation (5b) through the dependence of 
y on n c - The slight dependence of D &c on n c has already been neglected on deriving 
equation (5a) from equation (4). Instead of solving equation (9) and its coupled partner 
equation (12) by a direct numerical approach, an iteration process is used. This method 
allows y to be treated as a known "constant” at any iteration, since its value is available 
from the previous iteration. This procedure is not unlike the procedure used to treat 
nonlocal terms in the self-consistent field calculations in atomic physics. If y is 
"constant," then 


8 



f(r) = cosh yr + C 2 sinh yr 
Boundary condition (7b) relates Cj to C 2 as 
„ ^ yR tanh yR - 1 

L/ o — v_/ a ■ ■ 

i 1 tanh yR - yR 
Applying equation (7a) fixes -C^ as 


C l = 


U d x ' 

yR tanh yR - 1 , , 

cosh yr n + -r — r — n 5 — sinh yr n 

' 0 tanh yR - yR ' ° 


With Cj so defined, the steady-state solution may be expressed as 


n a" (r) * T-( cosh ^ + slnh >*) (10) 

Equation (10) simplifies whenever yR is either very large or very small. The magnitude 

of y is governed primarily by the local concentration of n„. To be specific, if yR £ 5 

yr ^ 

then tanh yR * 1, C 2 ~ -Cj, and Cj ~ r Q N a (T)e °. Then equation (10) is approximately 

I- 

n a °°(r) ~ ^ N a (T)e 7 ( °) (yR large) 


As a corollary large gradients in the distribution are to be expected whenever 
^a/ D ac R ~ 5 - 

To examine the limit as y — 0, rewrite equation (10) as 


na°°(r) 


r o^a^) (tanh yR - yR)cosh yr + (yR tanh yR - l)sinh yr 
r (tanh yR - yR)cosh yr Q + (yR tanh yR - l)sinh yr Q 


( 10 a) 


Then since for yR small sinh(yR) ~ yR + (yR)^/ 6 , cosh(yR) = 1 + (yR)2^2, and 
tanh(yR) ® yR - (yR)^/3 to order 3, equation (10a) becomes 
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«, g r 0 N a (T) -ry 3 R 3 / 3 + (y 2 R 2 - l)(yr - y 3 r 3 /s) 
r -r 3 R 3 /3 + (y 2 R 2 - 1 )(y? Q - 7 3 ^ 0 3 /^ 

_ r 0 N a (T) -r + y2( r 3/6 + rR 2 - R 3 /3) + ©(y 4 ) 

r -r 0 + y 2 ( r o 3 / 6 + r 0 R 2 - R 3 /3) + ©(y 4 ) 

then 


n a °°(r) =N a (T) (yR - 0) 

lim y— 0 

Expression for n c °°(r).- The and n £ solutions are coupled through the anni- 
hilation constant and may be solved by iteration. An explicit expression for the C con- 
centration is needed: 

D ac v2n c < ”( r ) ' Kc( n a) n c°°{ r ) = 0 < u ) 

By using the angle independent Laplacian and letting h(r) = rn c °°, equation (IT) can be 
reduced to 


- f|2|,(r) = 0 ■ (12) 

dr 2 

where 0 2 = K c ^jD ac . 

The boundary conditions analogous to equations (7a) and (7b) are 

n c °°(R) = N C (T) (13a) 

vS c “(r o ) = 0 (13b) 

where N C (T) is the number density estimated from measurements of outgassing pres- 
sure. Empirically, outgassing depends on time, temperature, and construction material. 
Equation (13b) implies that there is no net flux of C across the surface of the bead. 

With the same stipulations as applied to n a °°, the solution is 

h(r) = Cg cosh (3r + sinh /3r 
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where the boundary conditions require that 

tanh j3r 0 - 0r o 

vo — v/ 

0 * j3r Q tanh /3r Q - 1 


and 


C 4 = RN C (T) 


<3r Q tanh j3r Q - 1 


V v 

(tanh /3r 0 - ^r^cosh 0R + (j3r Q tanh /3r Q - l^sinh /3R 


Therefore, 


(tanh /3r Q - /3r 0 )cosh /3r + (/3r Q tanh /3r Q - l)sinh /3r 
(tanh /3r Q - |3r Q jcosh /3R + (j3r Q tanh /3r Q - l^sinh /3R 

Expressions for k & , and D ac .- Only binary diffusion is considered. Let 

represent the binary diffusion coefficient of species i and j. Then the lowest 
order approximation to the exact theory using hard sphere approximations for the colli- 
sion integrals gives (for example, see ref. 5, eqs. 7.4-38 and 8.2-5) 


n 


RN 

OO __ 


cn . = 3/7ikT\ 1/2 1 

J 8 \ m ij / pi + njjQ.. 


(15) 


where is the reduced mass and is the kinetic cross section for i - j 

collisions. 

Using perturbation expansions to solve the transport equations for multicomponent 
systems, the generalized diffusion coefficient may be written in terms of the binary 

counterparts. For example, if a third component k is present (eq. 11.2-41 of ref. 5), 
then 


D U 




n k( m k®ik/ m j ~ C ^ij) 

ni^Djk + nj^ik + “k^y 
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In particular, for the problem treated in this paper, 


^ac ^ac 


1 + 


n b( m b q) ab/ m c - ^ac) 
na^bc + n c CD ab + "b^ac 


(16) 


Note that when the concentration of n^, dominates, n a c Dj 3C and n c <r D a b are negligible 
in the denominator and 


D ac ~ ^ac 


1 + 




b u ab 
“c^ac 


- 1 = 


m b q) ab 


so that it is really the dominant species which hinders the flow of A as expected. 

As defined in the diffusion equations, the annihilation constant k is identified with 

the rate constant for a second-order chemical reaction. Kinetic theory considerations 

show k to be proportional to the cross section for encounter Q„„, the relative velocity 

1 /o 

of approach ocAr/m ac J ' , and the number density of the reacting partner. More detailed 
consideration introduces an activation energy AE representing an energy barrier which 
the reactants must surmount before a reaction can occur. Finally, an additional propor- 
tionality constant g accounts for the internal complexity of reactants in energy alloca- 
tion. In the appendix an expression for k is derived on the basis of Eyring's theory, 
where it is shown that 

K a,c = g Q ac( 8kT / 7rm ac.) 1/2n c,a c ' AE/kT (17) 

where A is atomic and 

g ~ 1 (C is atomic) 

0.01 = g = 0.1 (C is diatomic) 

0.0001 ^ g = 0.01 (C is polyatomic) 

Time -Dependent Solution 

Assume that conditions are such that equation (10) reduces to a decaying exponential 
at some intermediate radius r where r Q < r < R and that at time t = t Q , the source 
strength is instantaneously increased by a dimensionless factor A to the value 
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(1 + A)N a (T). Mathematically, this condition amounts to a source perturbation repre- 
sentable as AN a (T) 5^r - r Q j 6 (t - t Q ). In practice the surface of the alkali is dirty or 
contains an oxide coating so that the observed vapor pressure value is reduced to some 
fraction of the value (at temperature T) listed in the handbooks. If the surface is 
mechanically disturbed in some way, then some "clean metal" is momentarily exposed, 
and the vapor pressure is expected to jump discontinuously for a short time as described 
above . 

The time dependence is carried in a function ?j a (r,t) so that n^r,^ = n a °° + 77 a (r,t). 
Then for t > t , 

h a = "a 00 + = ( D ac v2 " “aH" + ( D ac v2 - « a K (r ’ t) 

In the preceding expression steady- state values of n a °° and (D ac v2 - k^i ^ 00 are zero. 
The function 77 must satisfy 

fj a (r,t) = (D ac v 2 - K a )j a (r,t) 

with 


77 a (r,t )=0 (tSt 0 ) 

Let rn = g and measure time from t_ by using time variable t = t - t . The func- 

d d U O 

tion g must satisfy 

cl 

d^g 

6 a (r ' T) ■ D ac -r - v a (r - r) 

dr 

This equation may be solved with Laplace transforms. The transform partners will be 
g a — - G a and r — s. Denoting the transform operator as then 

a {e a (r ' T >} - D ac £2 a (g a } - *a sf {sa (r ' T) } 

where the commutability of S£ and d2/dr2 has been used. From the properties of 
Laplace transforms (ref. 6 , p. 4) and the initial condition, 

sG a (r,s) - g a (r,0) = D ac G & (r,s) .- /c & G a (r,s) 
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where the steady-state value of g a (r,0) is zero. The solution is G a (r,s) = Cge^ r + Cge'£ r 
r 1/2 

with C = + s ]/^ac • The boundary conditions analogous to equations (7a) and (7b) 


are 


dGg 

dr 


(r,s) 


R ‘ G a (R ’ s) 


G a( r o’ s ) = r o AN a< T > 


and lead to 


-2?R 1 + R? 
^5 - - c 6 e TTr5 


r o aN a (T) 



1 a- Rt -g(2H-r 0 ) 
1 - RC 


Under the restricted range of applicability (£R large), 
c 6 “ ar o N a< T > e?r ° 

C 5 «0 

so that 

G a (r,s) ~ r Q AN a e ^ ^ 

The inverse transform may be obtained by using property 3 and transform number 85, 
respectively, on pages 243 and 250 of reference 6 as 


ST^Ga^s)} = r o AN a (T) 


>a T ( r - r o) e 


■(r- r o) 2 /47 D ac 


K^ac) 


\l/2 
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In terms of original variables 



Ar o( r - r o) N a( T ) e 


- K a( t - t o)-( r - r o) 2 / 4D ac( t - t o) 





?7 a (r,t) = 0 



It should be possible under certain conditions to observe responses to changes in 
the steady-state distribution, the time dependence of which is carried primarily by 
exp(-/c a T). Time-dependent observations may help to determine at least the magnitude 
of x a when its direct calculation is not feasible. 

Numerical Example of Steady-State Solution 

A numerical example will help to describe the method and point out the 
salient features of the solution. If' only roughing pumps are used in preparing 
the absorption vessels, residual total pressures may be as much as 10 pm Hg 
(1 (im Hg = 10-3 torr). When alkali is introduced, it in effect acts as a ’’getter" 
for many of the oxygen- containing outgassing constituents. This process removes 
atomic phase alkali. Outgassing rates depend markedly on temperature and time 
history of heating, the latter because there is essentially only a finite amount of 
material to outgas. For the purpose of identifying the level of contamination in this 
example with a single parameter (time), it is assumed that p = 10~2 exp(-t/q) where p 
is in torr, t is in hours (1 hour = 3600 sec), and q acts as the pressure time con- 
stant. This relationship oversimplifies the actual behavior and overestimates the con- 
centration of harmful constituents, for example, oxygenated species. If q is 3 hours, 
then this value is consistent with data in figure 4.5 of reference 7 for total outgassing 
and desorption of nonpolished aluminum chambers at 200° C (472 K) after Ng flushing. 
The relationship is used for times greater than or equal to 2 hours. This long time con- 
stant would at least be doubled at room temperature on the basis of results in reference 7. 
Typical diatomic molecular diameters are 3A to 4A (lA = 10“ ^ m ). For calculating the 
cross sections in this example, 3.5A is used as a mean diameter and the mass of C will 
be taken as 30 amu. 

Species A will be Cs and B is taken to be Xe at 400 torr. Vapor pressure 
data for Cs is taken from page 2425 of reference 8, and the atomic and molecular radii 
for cross-sectional calculations from reference 9, pages 45 to 47. The activation energy 
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is assumed to be zero since Cs reacts violently with oxygenated compounds even at 
room temperature. Table I assembles the input data for the example. 

Equations (10) and (14) are solved by iteration. An initial guess for n °°.(r) is 
made to allow K a to be estimated from equation (17). The expression for y is 
calculable and hence n °°(r) to first approximation is obtained from equation (10). The 

a. 

quantity k and hence /3 may be used in equation (14) for a new estimate of n c °°(r) 
and the procedure repeated until selfconsistency is achieved. Convergence usually occurs 
within 10 or 20 iterations and, as expected, is slowest in the regions of strong gradients. 

TABLE I.- VALUES USED FOR NUMERICAL EXAMPLE 

2.20 x 10- 22 g 

2.18 x 10-22 g 

. . . . 5.00 X 10-23 g 

473 K 

0.5 cm 

10.0 cm 

. 1.42 x 10 19 cm -3 

2.48 x 10 15 cm' 3 

3.53 x 10 14 exp(-t/10 800); t in s 

83.0 x 10- 16 cm2 

91.0 x 10-16 cm 2 

44.0 x 10' 16 cm2 

0.0 

0 . 01 , 0.001 

Figure 3 shows computed profiles as the concentration of contaminant depletes in 
accordance with the simplified time dependence for outgassing and desorption listed in 
table I. The time intervals are also comparable to run times encountered in the absorp- 
tion experiments under consideration. The figure is only intended to represent the sharp 
demarcations that result from this approach. Exact values for the dependence at the reac- 
tion zone would require identification of species C and calculation of the associated par- 
tition functions, and so forth. The demarcation becomes less abrupt as g becomes 
smaller, which is to say as C becomes more complex. It then takes more encounters 
to insure a reaction. The example seems to indicate that spatial gradients appear when- 
ever N c (T) levels approach within three orders of magnitude of N a (T). Use of diffu- 
sion pumps to prepare the experiments would result in approximately four orders of mag- 
nitude reductions in the initial N C (T) levels at the outset and should adequately remove 
the gradients in small cells loaded by distillation techniques. For large cell loading tech- 


m b 

m c 
T . 


R. . 

n b • 
N a (T) 

N C (T) 

Qac 

Q ab 

Q bc 

AE 

g 
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niques, as used in reference 1, the alkali may be used as a viable getter if the alkali sur- 
face can be maintained in a "clean" state. 

Figure 4 shows the effect of reducing the buffer gas pressure by a factor of 50 and 
is interpreted to mean that the A atom has less interference as it moves away from r Q 
and less chance of colliding with C so that the demarcation zone moves toward the wall. 
The experimental evidence (ref. 1) also seems to indicate this result. 

CONCLUDING REMARKS 

By working with a simplified model it has been shown that the mutual diffusion of 
chemically reacting species in a medium of inert gas atoms can result in severe spatial 
gradients of both reacting species. If the alkali annihilation constant is the second- 
order rate constant for the reaction of A with C , then rather severe gradients in the 
concentration of A are expected if ^ a ^/D ac R * 5 where D ac is the multicomponent 
diffusion coefficient and R is a characteristic dimension separating the sources of the 
reactants . 

A numerical example for a 400-torr bath pressure showed that when the source 
strength of one of the reactants approached to within three orders of magnitude of the 
other, severe spatial gradients of both reactants resulted. This result is germane to 
some optical absorption measurements with alkali metals, since optical absorption tubes 
frequently operate in a temperature range that puts the alkali vapor pressure within three 
orders of magnitude of outgassing levels for oxygenated compounds from the vessel wall. 
Recent experiments are cited which, in fact, support the predictions of the present model. 

The analysis also indicates that the gradients become less severe as the complexity 
of the reactants increases; for example, gradients are more severe for atom-atom cases 
than for atom -polyatom situations. 

The effect of the inert gas bath pressure is carried through the use of a generalized 
diffusion constant. Decreasing the bath pressure makes the reaction zone less abrupt and 
moves it toward the boundary supplying the minor reactant - a result also compatible 
with observation. 

Perturbations in the source strength are shown to propagate toward the reaction 
zone and to relax there with a time constant inversely proportional to the rate constant 
for the reaction. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 5, 1974. 
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APPENDIX 


ESTIMATION OF k 0 

d 


The experiments usually neither measure nor identify species C. A reasonable 
estimate may be given, however, based only on knowing whether the species is atomic or 
polyatomic. Likely outgassing candidates, when aluminum or stainless steel vessels are 
involved, might be O, C^, HgO, CO 2 , CO, and some hydrocarbons from O-rings. 

Assume for simplicity that every encounter of A and C removes A as an 
atomic constituent. Then the loss rate of n a can be estimated with simple kinetic theory 
arguments. For a Maxwell- Boltzmann gas, the mean relative speed of A and C is 
given by ^kT/irmac^/^ where m ac is the reduced mass of the pair. An A atom 
will, on the average, sweep out a volume Q ac ^8kT /7rm ac j*/^ per unit time where Q ac 
is the cross section for encounter. Multiplying by the number densities of A and C 
gives their mean collision rate and a crude estimate for the loss of 



(18) 


In chemical kinetics this condition would be called a second-order reaction because 
it is of first power in each concentration. The term in brackets is called the rate con- 
stant. A generalization of the rate constant will be based on chemical kinetics. (See 
chapter 3 of ref. 10.) Expression (18) will be modifed by a fraction g which is a func- 
tion of the complexity of C and by an additional temperature-dependent term. 

Much of reaction rate theory is based on the assumption that the colliding partners 
comprise an activated complex which has all the features of a stable molecule save one - 
that one of its vibrational modes dissociate the complex into its final products. In fact, 
the frequency of this special vibrational mode is identified with the rate of reaction. 
Furthermore, both the initial and final products are assumed to be independently in 
equilibrium with the activated complex. The quantity X* is the concentration of acti- 
vated complex. Mathematically, 

n a + n c Z X* Z Final products 


Let Z represent the partition function. Statistical mechanical arguments show that the 
relative concentrations of various species are related by 


X* _ Z* -AE/kT 
n a n c z a z c 


( 19 ) 
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where AE is the "activation energy” and is defined as the average energy of the mole- 
cules which react minus the average energy of all the molecules. (For example, see 
Rice, ref. 11, p. 513.) Pictorially, AE represents a kind of "energy hump," which the 
colliding partners must overcome to be considered a complex, and empirically AE is 
responsible for the strong temperature dependence of most rate constants. 

If the frequency of the vibrational mode which leads to dissociation is r( d ), then the 

contribution of this mode to Z* will be |l - exp^-hr(^)/kTjj . Since z/( d ) belongs to 
a very weak coupling, it will be very small with respect to kT/h which is approximately 
6 x 10 12 sec - l at 300 K. Consequently, the exponential is representable by the first two 
terms in its Maclaurin expansion. The partition function to this accuracy is kT/hi'W. 
Designate the partition function with this factor removed by Z*; that is, 

i 

Z* = Z^hv^/kr) (20) 

Substitute equation (20) into equation (19) and rearrange to obtain 


xM d ) 


kT Z * e -AE/kT 
h z a z c 


n a n c 


( 21 ) 


which will reduce to equation (18) if C is also an atom. 

To the extent that energy in a molecule can be identified as translational, rotational,, 
vibrational, or electronic, it is possible to factor the partition function. Let the partition 
function for a single mode be z with the superscript identifying the type; then, 


w = (27rmkT) 1//2 / h 

(22a) 

W = (8ir2ikT ) 1 / 2 /h 

(22b) 

( y ) = [l - exp(-ht'/kT)] -1 

(22c) 

( e ) = Degeneracy of lowest state 

(2 2d) 


where I in equation (22b) is the moment of inertia, equation (22c) is the harmonic oscilla- 
tor approximation, and equation (22d) is true only for temperatures below 1000 K. Order 
of magnitude estimates (in cgs units) for the individual modes are 10^ for z(t), 10 to 
100 for z( r ) and 1 for z(v) and z( e ). It is convenient to define z( v )/z( r ) = e; then 
0.01 ^ e ^ 0.1. 


1914 



APPENDIX 


If A and C are atoms, X* is a diatomic molecule. The electronic mode is 
neglected. Each atom has three degrees of translational freedom. The quantity Z* has 
its only degree of vibrational freedom removed to form Z *; therefore 


Z a = z a ^ 3 ~ z^ 3 


Z c = z c (t)3 «z (t)3 


Z* = z W 3 z W 2 *z(t) 3 zW 2 


Substituting these relations into equation (21) 


r 


X V d > = ) 




kT z * 


(t)3 (r)2 


h (t)3 (t)3 

z a z c 


,-AE/kT 


n r) n ; 


(23) 


If the z terms are calculated with equations (22) and it is noted that I = Q ac 

then equation (23) reduces exactly to equation (18) except for the activation energy term. 

Identify g times the expression in braces with x a . Then g is 1 if C is atomic. 

If C is diatomic, then X* is triatomic less one degree of freedom 


and 




Z. . 


Z. =z< t ) 3 z( r ) 3 z< v ) 2 


Z,Z 


a-c 


atom -diatom 


lZ a Z _/ 

' '''atom -atom 
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If C has three or more atoms, g is e^. In summary 

“a = ^ac( 8kT / 7rm ac) 1/2n c e " AE/kT 

where g~l if C is atomic; O.Olgggo.l if C is diatomic; and 0.0001 = g = 0.01 
if C has more than two atoms. 
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ISOTHERMAL ENCLOSURE 


^£1 


ABSORPTION TUBE 






7 


Cs RESERVOIR 


OPTIC AXIS 


(a) Schematic of a typical absorption cell set up to indicate the location 
of Cs reservoirs with respect to optical axis. 



POSITION ALONG OPTIC AXIS , cm 


(b) Data showing the Cs vapor distribution above two separated reservoirs 
for loading pressures of 10, 100, and 500 torr He. Measurements are 
made along the optical axis depicted in figure 1(a). 


Figure 1.- Experimental alkali concentration measurements from reference 1. 






RADIUS, cm 


RADIUS, cm 


Figure 3.- Steady state distributions (a) to (d) which result when the contaminant level 
n c changes at the boundary r = R while the alkali density n a is held at the 
200° C vapor pressure value at the bead's surface r = r Q . Contaminant level is 
parametrized by time according to the simplified model in table I. Buffer gas 
pressure is 400 torr. 
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